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production processes is discussed and confronted with experimental data on e + e - annihilation processes, pp and pp scatterings and 
heavy ion collisions. 

PACS: 02.50.-r; 05.90.+m; 24.60.-k 

Keywords: Information theory; Nonextensive statistics; Thermal models 



> 

m 

(N 

m 
o 

a: 
a: 

1 Introduction 



High energy multiparticle production processes are most naturally described by 
statistical models (see [U[2j for a historical background and 011! for the most re- 
cent developments). Their results are usually interpreted in the thermodynami- 
cal sense, with temperature T and chemical potential fi entering with their usual 
meaning. However, it has been recognized that in this branch of physics one very 
frequently encounters "thermal-like" form of distributions in some variable x, say 
oc exp(— x/T), without system under consideration being in any kind of thermal 
equilibrium pj. It is enough that out of the huge number of produced secondaries 
only some part is registered by detectors, and out of them only one or two are 
selected for final scrutiny. The averaging emerging this way is then equivalent to 
the action of some "heat bath" characterized by parameter T. Once this is real- 
ized it is then obvious that there are situations in which such " heat bath" is more 
complicated (for example nonextensive) and needs additional parameter (s) to be 
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described properly PEUH]. In this way the nonextensive Tsallis statistics charac- 
terised by the non-extensivity parameter q enters in a natural way Pll0pillf l2j. 
On the other hand, such apparently "thermal- like" behaviour of some distribu- 
tions arises very naturally in many physical applications of information theory 
(MaxEnt) [T3] (with Shannon form of the corresponding information entropy, its 
nonextensive form uses Tsallis entropy characterised by the same parameter q as 
mentioned above, such that for q — > 1 one recovers the usual Shannon entropy). 



2 Results 



The usefulness of information theory is most obvious in situations when one has 
to "guess" the most probable (and least biased) distribution p(x) of some quantity 
x using only limited amount of information given in terms of finite number n of 
some observables: Rk=i,..., n = {Rk{x)) {p{x) is normalized to unity and average 
(. . .) is with respect to p(x) in extensive and to [p(x)] q in nonextensive approaches, 
respectively [T3|I9] . In high energy multiparticle production processes it allows to 
find, in a maximally model independent way, the real information content of the 
experimental data considered |TH[TK1[TH1[TTI 18]. From (Hj we know therefore that: (i) 
- particles produced in high energy collisions are mostly located in one dimensional 
phase space, i.e., they have limited transverse (with respect to collision axis) mo- 
menta px {{pt) is finite) and are fully characterized by their longitudinal momenta 
Pl = Ursinhy ; (ii) - only fraction K G (0, 1) (called inelasticity) of the initially 
available energy W is converted into produced particles. As was shown recently by 
us [IB] the fact that multiplicity distribution of produced secondaries is no longer 
of Poissonian type but follows much broader Negative Binomial (NB) form can be 
accounted for only by using nonextensive approach with q given by the parameter 
k of NB, q = 1 + which represents the amount of dynamical fluctuations in 
the number of produced secondaries [18J. 

We are usually interested in single particle rapidity distributions, which in the 
information theory approach are given by 

1 dN 1 

P(v) = -^^T = -7r ex Pq(-Pq ■ Vrcoshy) . (1) 



N dy Z t 



* Here ht = \J li 1 + (pt)' 2 is the so called transverse mass of particle with mass \x and transverse momen- 
tum pt and y denotes rapidity of the particle, variable defined in such way that its energy is E = cosh y. 
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This form is identical with that used in statistical models but now Z q and (3 q are no 
longer free parameters to be fitted when comparing with experimental data but in- 
stead are given by the normalization condition and energy conservation constraint, 

Y T Y P k ■ W 

J dyp(y) = 1 and J dy /i T ■ cosh y ■ [p(y)] q = q (2) 

— V —Y 



(where ±Y m are maximal rapidities available in rest frame of hadronizing source, 
see P^UTH] for details). There are therefore two parameters: nonextensivity q re- 
sponsible for dynamical fluctuations and (/-inelasticity K q to be deduced directly 
from data essentially in a model independent way . This method works perfectly 
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Fig. 1. Examples of applying eq. Q to: (a) rapidity spectra for charged pions produced in pp and pp 
collisions at different energies |19|2()l21l22j : (6) similar data obtained for the most central Au+Au collisions 
[2"3] : (c) rapidity spectra measured in e + e~ annihilations at 91.2 GeV [H] (dotted line is for K q = 1 and 
q = 1 whereas full line is our fit with K q = 1 and q = 0.6). 

well for pp and pp collisions [T8| . cf. Fig. [T^l (where it provides us with the first 
model independent estimations of energy dependence of the mean inelasticity pa- 
rameter and with its distribution). Here we show that it works also quite well for 
similar data on rapidity distributions obtained in Au + Au collisions, cf. Fig. [Hd. 
The Au + Au data are for the most central events (covering collisions proceeding 



Here exp g (x/A) = [1 + (1 - q)x/A] 1/{1 - q) q -- 
ing of K q for q ^ 1 case is solved by noticing 
physical meaning, is K q = K q j (3 — 2q). 



exp(x/A)). Possible question concerning physical mean- 
that the inelasticity parameter in this case, which has 
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with impact parameter range — 6%). They can be fitted choosing K q = 1 and 
then q = 1.29, 1.26 and 1.27 for energies 19.6, 130 and 200 GeV, respectively 
(the (/-inelasticity was therefore equal to K q = 0.42, 0.48 and 0.46) ** * . The most 
interesting are, however, results for e + e _ annihilations (cf. Fig. Hh) for which, by 
definition, K q = 1 (i.e., always all energy of initial leptons is available for the 
production of secondaries) and which can be fitted only with q < 1 (in our case 
q = 0.6). 

This point deserves closer scrutiny. One could argue that because fit in Fig. [Tb is not 
perfect (there are some discrepances for small rapidities and there is a tail at large 
values of y) there is nothing to be said before they are not addressed. But results for 
q = 1 clearly show that these discrepancies are not connected with the particular 
value of q but rather with some additional mechanisms operating here action of 
which would, however, change our results only slightly (for example, a possibility 
of two rather than one source or ^/-dependent (pr), as mentioned already in 
With the above reservations let us then take a closer look at the possible origin of 
q < 1. We have already encountered similar situation when in (T7j we have fitted 
single particle distributions assuming implicitely that K q = 1 and discovering then 
that one could get fairly good agreement with data for q < 1 only. That was because 
in this case only q < 1 leads effectively to cutting-off a part of the phase space 
(once it is taken too big) mimicking therefore action of the inelasticity parameter 
(cf. |T8|). On the other hand, we know that q ^ 1 signals presence of fluctuations 
in the system PTITTUT^ and is given by normalized variations of these fluctuations, 
in our case they would be fluctuations of temperature T = 1/(3 parameter. So far 
it was widely discussed only for the q > 1 case fTT|[T2] but formally it covers the 
q < 1 case as well, cf. |Hj. However, in this case temperature T does not reach an 
equilibrium state because in this case 

T = T - (1 - q)E (3) 

instead reamining constant, T = To, as is the case for q > 1. In this case we have 
a kind of dissipative transfer of energy from the region where (due to fluctuations) 
the temperature T is higher (for example, in our case from the quark (q) and 
antiquark (q) jets formed in the first e + e _ — > q + q to gluons and qq pairs and later 

*** Actually data in Fig. are presented for the so called pseudorapidity r\ defind not by energy E 
and longitudinal momentum pi as is the case of rapidity y but by total momentum p and longitudinal 
momentum p^ instead. 



4 



on to finally observed hadrons). It means therefore that q < 1 signals that in the 
reaction considered, where K q = 1 and we have to account for the whole energy 
exactly, conservation laws start to be important and it is not possible for stationary 
state with constant final temperature to develop but temperature T depends on 
the energy * * ** and for large energies tends to zero (notice that from eq. (JHJ) one 
has limitation on the allowed energy of the secondaries: E <Tq/(1 — q))** . 



3 Summary 



We have demonstrated that information theory (know also as MaxEnt method) can 
play very important role in high energy physics providing simple, highly model 
independent, estimations of single particle distributions and allowing to reliably 
estimate the amount of information provided by experimental data. Its nonexten- 
sive version extends applicablity range of MaxEnt by including also some intrinsic 
fluctuations in the hadronizing system visible as broadening of the multiplicity dis- 
tributions. As was shown here, single particle distribution data on all kinds of col- 
lisions, starting form very elementary e + e _ annihilations, via pp and pp collisions, 
and ending with the most complicated heavy ion scatterings, can be described by 
formula (CD) depending on very limited number of parameters: inelasticity K and 
nonextensivity parameter q. Fairly good fits were obtained in all cases. To proceed 
further one should concentrate now on these parts of the phase space where dis- 
crepancies occur and then introduce one-by-one some additional hypothesis and 
check whether they lead to the better agreement with data (see, for, example |25|). 
Such approach allows to avoid assumptions, which whereas looking promising, are 
not justified. 
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Actually, when fluctuations depend on it in the same way the relative variance u remains constant 
and q = 1 — u, cf. [TT] , 

★ * * * * Notice that analysis of distribution in the same process using a kind of q- version of Hagedorn 
model |26j reports q > 1 instead. This is, however, what we call qx in and this is different from q 
considered here (called q^ in where we compared both q's for pp collisions). The reason for such 
different behaviour is that pt's considered there are essentially not influenced by conservation laws but, 
in our language here, reflect instead a kind of stationary state with q > 1 and energy independent T. 
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